Quantum and classical localization in the lowest Landau level 



cn 
o 
o 

(N 

5h . 

< 
OO 



X3 

o 



> 

o 
m 
o 

-I— > 



o 



X 
S3 



Nancy SandlerQ Hamid R. Maei, and Jane Kondev 
Department of Physics, Brandeis University ,Waltham, MA 02454 
(Dated: February 2, 2008) 

Spatial correlations of occupation probabilities, if their decay is not too fast, can change the 
critical exponents for classical percolation. From numerical studies of electron dynamics in the 
lowest Landau level (LLL) we demonstrate the quantum analog of this effect. Similar to classical 
percolation, we find that the extended Harris criterion applies to localization in the LLL. These 
results suggest experiments that might probe new quantum critical points in the integer quantum 
Hall setting. 



Quantum or classical dynamics can lead to localization 
of a particle moving in a random environment. The lo- 
calized state in the quantum case is described by a wave- 
function which spreads over a finite distance, the local- 
ization length. In the classical case particle trajectories 
are closed orbits limited in extent. The mechanism by 
which localization can occur in the two settings is very 
different as quantum interference, which plays a crucial 
role in quantum localization, has no classical counter- 
part. In this letter we report on the intriguing relation 
between classical and quantum localization in the setting 
of the integer quantum Hall (IQH) effect, when a power- 
law correlated disorder potential is present. 

The problem of localization in the IQH setting has been 
studied since the effect was first discovered. Localization 
is believed to be responsible for the plateaus in the Hall 
resistance which are observed as a function of the applied 
magnetic field 1]. Transitions between adjacent plateaus 
have been identified with a zero-temperature quantum 
critical point 2] which is characterized by the divergence 
of the localization length as the magnetic field is tuned 
to its critical value. Support for this picture is provided 
by observations of critical scaling in experiments of Wei 
et al.^ and Koch et al.\^, which found a localization 
length exponent Vq « 2.3. 

Theoretical studies of localization in the IQH system 
usually take the semiclassical picture ^] as their starting 
point. Here the electron's motion is split into the fast 
cyclotron rotation and the slow E x B drift of the guiding 
center along lines of constant electrostatic potential. The 
potential V^(r), caused by impurities, is assumed to be 
random. Quantum mechanics enters at the saddle points 
of the potential where the electron can tunnel from one 
equipotential orbit to another, while the perpendicular 
magnetic field results in a random Aharonov-Bohm phase 
between two tunnelling events. A lattice model which 
describes the network of saddle points was introduced by 
Chalker and Coddington and it leads to a value of 
the localization exponent consistent with experiments. 

The classical limit of the network model corresponds 
to percolation 0. The classical trajectories are equipo- 
tential lines of ^(r), which can be mapped to the hulls of 
percolation clusters However, little is known about 
quantum localization when the classical limit is described 



by correlated percolation. Namely, as shown by Wein- 
rib 9j, power-law correlations of the random potential 
can change the percolation critical point leading to a dif- 
ferent fractal geometry of its equipotential lines. This 
happens when the exponent a, which characterizes the 
spatial decay of the potential correlations (see Eq. EJ, 
satisfies a < 2/^^, i^c = 4/3 is the correlation length 
exponent for percolation [lOj. Here we investigate the 
effect that changing the fractal geometry of classical tra- 
jectories, has on quantum localization in the IQH system. 
We find new localization critical points when a < 2/i/g, 
i.e., the quantum analog of correlated percolation. (This 
problem was previously studied by Cain et al. using real 
space renormalization group techniques, but no definite 
proof of the effect was found llj.) 

Quantum dynamics in the lowest Landau level The 
Hamiltonian for two-dimensional spinless electrons con- 
fined to the x-y plane, and under the combined effects of 
a magnetic field B = Bz and a random potential V{y), is 
H = I/277i(p - eA)^ -I- y_k Pk- Here A is the vector 
potential, pk = e**' '" is the one-particle density operator, 
and 14; is the Fourier transform of the disorder poten- 
tial. At high magnetic fields (or low temperatures), the 
quantum dynamics of the electron are governed by the 
projection of the Hamiltonian onto the lowest Landau 
level (LLL), 
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where p is the projected density operator. We focus 
our attention on the localization properties of the eigen- 
f unctions of H . Previous numerical studies indicate 
that the localization length increases towards the cen- 
ter of the LLL band, with a localization length exponent 
Uq K, 7/3. In these studies the random potential was as- 
sumed to be short range correlated. Here we take up the 
question of how power-law correlated potentials might 
affect the value of Vq. 

Recently, Sinova et al. 13] have shown that the local- 
ization length exponent can be computed from the dis- 
order averaged density-density correlation function pro- 
jected onto the LLL. It was further shown by Gurarie and 
Zee ^3 that the classical limit of the time evolution of 
this correlation function describes electron drift along the 
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equipotentials of y(r), as expected. These results lead to 
a simple scaling argument for the subdiffusive spreading 
of wavepackets in the lowest Landau level . Namely, 
start by assuming that the disorder potential completely 
breaks the degeneracy of the LLL, leaving only one ex- 
tended state at energy Ec — 0. Let us examine the time 
evolution of a wavepacket {ipEit)) constructed from lo- 
calized eigenstates of Eq. Q that are taken from an in- 
terval of width A centered around E. The dispersion of 
this wavepacket, < Ax'^{t) >e= {ipEit)\x'^\ipE{t)), is 
expected to be diffusive for short times {t <^ (,'^{E)/D, 
D— diffusion constant), attaining a constant value set by 
the localization length £,{E) at long times. This can be 
summarized by a scaling form: 

<Ax^t)>E = Dtfl^-^^ (2) 

where the bar denotes disorder averaging and / is a scal- 
ing function, with properties: f{x) — > const for x <^ 1, 
and f{x) const/a: for a; ^ 1. Now, if we construct 
a wave-function in the LLL as a superposition of the 
wavepackets \^E{t)) its dispersion can be written as 

< Ax^t) > ~ < Ax2(t) >E, (3) 

where the contributions of the non-diagonal terms can 
be neglected since only the state at = is extended. 
Replacing < Aa;^ (t) >£ by Eq. Q , making use of ^ {E) ~ 
E^"^"! , and the non-singular nature of the density of states 
for i? ^ 0, we obtain: 

< Aa;2(t) > - {D i)i-i/2^, (4) 

Therefore, the spread of the wavefunction at late times 
is subdiffusive due to the presence of the extended eigen- 
function at the band center, and the anomalous diffusion 
exponent 9 = 1 — l/2vq. This forms the basis of the 
numerical method we use below to compute Vq. 

Power-law correlated disorder potential To study 
quantum dynamics and localization in the LLL in the 
presence of power-law correlated disorder potentials, we 
make use of the Hamiltonian in Eq. with V{k) a ran- 
dom Gaussian variable with zero mean and variance given 
by |k|"~2^ After Fourier transforming this leads to a po- 
tential which is power-law correlated in real space: 

{V{v)vm^^. (5) 

The projected Hamiltonian matrix is written in the basis 
of states given by the Landau eigenstates on the torus, 
which in the Landau gauge can be written as: 

oc 

'4}i{x,y) = ^ exp(27ri|-(iVm + 0) 

m— — cxj 

xexpf-^^[x-4(/ + 7Vm)]2) . (6) 



Here I goes from to A'' — 1, labelling the N states on the 
torus, N is the number of flux quanta through the torus, 
and L is the system size. These wavefunctions are defined 
on [0, L) X [0, L), they are centered at a; = Ll/N and have 
appreciable amplitude in a narrow strip of width given 
by the magnetic length tc. = \/ L"^ /'inN . (Below we set 
i = 1, rendering {x,y) dimensionless.) 

The electron density operator p{ki,k2) = 
exp(27ri(fcia; + ^2?/)) projected onto the LLL is a 
matrix of the form: 

p{ki,k2) = expf ^^-^TT j L(fci,fc2), (7) 

where fci, fc2 are integers in the units chosen, ((/ci, k2) ^ 
(0, 0)), and the matrix elements of L{ki, ^02) are given by: 

[L{k,,k2)],^ = e2-['=^ k./2+k,iH~WN |_^^ 

(8) 

(The formalism used to project the density operator pk 
onto the LLL was developed in Ref. HH ) 

For a given realization of V{r), the eigenvectors and 
eigenvalues of H are obtained by exact diagonalization 
|l7l |. A wave-packet localized along x is made using all 
the eigenstates. As it evolves it spreads, and its width 
in the x-direction is computed as a function of time. For 
one realization of the short-range correlated random po- 
tential, Fig^ shows a set of snapshots for the evolution 
of a wavepacket initially localized at a: = 0.5. To ob- 
tain < Ax^ (t) > the width of the wavepacket is averaged 
over all initial positions and over 1000 different disorder 
realizations; Eq. ^is then used for computing Vq. 

First, we checked this numerical procedure for a short- 
range correlated disorder potential, where the Fourier 
components of V{r) (in units of huc) are independent 
and uniformly distributed on the interval [—0.5, 0.5]. The 
value of 61 = 1 — 1 /2i'q is given by the slope of < Aa;^ (t) > 
when plotted on a log-log graph. To take into account 
finite size effects we repeat the computation of 9 for 
200 < A^ < 1500 and plot it as a function of 1/A^; see in- 
set of Fig. [3 By linear extrapolation we find i^g = 2.33(9) 
in the A^ — + 00 limit, consistent with previous results. We 
then repeat the computation for long-range correlated 
disorder potentials for different values of the parameter 
a. Results are shown in Fig. |5| 

For each < Ax'^{t) > curve in Fig. 12 three regions 
can be identified. Due to the finite system size, at long 
enough times the spread of the wave-packet reaches a 
constant value which corresponds roughly to 70% of the 
system size. The critical region corresponds to interme- 
diate times, while the short time behavior is described 
by the slope 9 k2. Ballistic motion {9 = 2) follows from 
aperturbative calculation of < Ax'^{t) > at short times 
|l8| . Comparison among different < Ax'^{t) > curves 
shows that the value of the slope in the critical region 
begins to increase for a < a* ~ 0.75 — 0.80. 
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FIG. 1: Evolution of a wavepacket composed of all eigenstates and initially localized ai x = 0.5. The real function plotted is 
the square of the amplitude of the wavepacket. 
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FIG. 2; Disorder averaged spread of the wave-packet, for 
A'^ = 1000 basis states in the LLL, as a function of time, 
and for different disorder potentials. Anomalous diffusion is 
evident at intermediate times. Inset: The exponent 6 for 
different values of A'^ in the case of a short-range correlated 
disorder potential. 



Extended Harris criterion The effect of short-range 
correlated disorder on a critical point is summarized by 
the Harris criterion The criterion states that criti- 
cal exponents for the disordered and the clean system re- 
main equal as long as the value of the correlation length 
exponent v satisfies; dv ~ 2 > 0, where d is the dimen- 
sionality of the system under consideration. It is derived 
by demanding that the fluctuations of the random po- 
tential within a volume set by the correlation length do 
not grow faster than its average value, as the transition 
is approached. An extension of the criterion was pro- 
posed by Weinrib and Halperin to include power- 
law correlated disorder potentials, like that in Eq. ||SJ|. 
When applied to two-dimensional percolation this "ex- 



tended Harris criterion" reads [9j: 

2 2 

a > a* ^ Vc = — ; a < a* ^ — — ■ (9) 
a* a 

where — 4/3 and a* = 3/2. However, there is to our 
knowledge no evidence of the validity of the extended 
Harris criterion for quantum critical points. Since the 
localization transition in the LLL is believed to be the 
quantum counterpart of the classical percolation transi- 
tion, it provides an ideal testing ground for the applica- 
bility of the criterion in the quantum realm. 

With this in mind we calculate the value of the anoma- 
lous diffusion exponent 6 — 1 — l/2vq by computing the 
value of 9 in the critical region for different values for 
the degeneracy of the LLL, ranging from N = 330 to 
N — 1000, and then extrapolating to cxd. Results 

from this computation, as well as the theoretical predic- 
tion based on the extended Harris criterion are shown in 
Fig. |31 In the quantum case the criterion reads as writ- 
ten in Eq. ® except now a* = 2/z/, = 0.786(8). In the 
inset of Fig. Owe also show results of numerical analysis 
of classical electron motion in the same long-range cor- 
related random potential, and com pare to the prediction 
from the extended Harris criterion The classical re- 
sults were reported previously in Ref. [2l| . Our data con- 
firm the validity of the extended Harris criterion for both 
quantum and classical electron dynamics in the LLL. 

The main result of this paper is contained in Fig. 13 
From the striking resemblance between classical and 
quantum behavior, we conclude that disorder correla- 
tions affect in a qualitatively similar way quantum and 
classical percolation. Indeed, the quantum version of the 
extended Harris criterion can be argued in close anal- 
ogy with its classical counterpart. In this case we should 
consider the fluctuations of the random potential in a 
volume set by the localization length £,{E) ^ £^"1 , and 
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FIG. 3: The anomalous diffusion exponent 9 — 1 ~ l/2uq 
characterizes the spread of a wave-packet in the LLL. The 
power a describes the long-range nature of the correlations 
of the random potential. Full lines are theoretical curves ob- 
tained from the extended Harris criterion. The inset shows 
analogous results for the case of classical electron motion [l^ . 

compare them to E. The calculation that follows llS l 
exactly parallels that done by Weinrib and Halperin [2flj , 
and Eq. Q follows with Vq replacing Vc- 

An interesting consequence of the applicability of the 
extended Harris criterion to the IQHT, is the conclusion 
that the simple relation between the classical and quan- 
tum localization length exponent, i^q = i^c + ^ |22l |. can 
not be valid in general. According to the quantum ver- 
sion of Eq. , for a long-range correlated disorder with 
0.8 ^ a < 1.5 the quantum localization exponent will be 
the same as in the case of short-range correlated disorder 
while the classical one will vary continuously with a. 

We conclude by speculating about a possible experi- 
mental test of our findings. The random potential present 
in samples that exhibit the IQH effect is thought to be 
short-range correlated, due to the nature of the spatial 
distribution of impurities within the heterojunction 0. 
However, the effect of a random magnetic field in the 
presence of a much stronger constant rnagnetic field is 
equivalent to that of a random potential 23l| . This opens 
up the possibility of engineering the properties of the dis- 
ordered environment seen by the electrons, by applying 
a random magnetic field. This was, for example, demon- 
strated in Ref. ,24 using a magnetic material with a rough 
contact surface placed in close proximity to the electron 
layer. The height fluctuations of the material surface, 
and hence the correlations in the random field, can in 
principle be controlled by adjusting the growth condi- 
tions under which the material is made. Magnetic dec- 
oration techniques or magnetic force microscopy can be 
used to determine spatial correlations of the field and 



thus extract the value of a. At T ~ 25mK, the satu- 
rated width of the Hall conductance step is AB ~ 0.5T 
for AlGaAs/GaAs structures. The temperature and 
AB place limits on the size of the random field fluctu- 
ations, which should roughly be in the interval lO^'^T 
< SB < lO^^T . Whether or not magnetic materials 
with these properties can be prepared, and placed in 
close proximity to the two-dimensional electron gas so 
as to effect its dynamics, remains to be seen. 
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